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ABSTRACT 

NGC 1313 is the most massive disk galaxy showing a flat radial abun- 
dance distribution in its interstellar gas, a behavior generally observed 
in magellanic and irregular galaxies. We have attempted to reproduce 
this flat abundance distribution using a multiphase chemical evolution 
model, which has been previously used successfully to depict other spi- 
ral galaxies along the Hubble morphological sequence. We found that 
it is not possible to reproduce the flat radial abundance distribution in 
NGC 1313, and at the same time, be consistent with observed radial dis- 
tributions of other key parameters such as the surface gas density and 
star formation profiles. We conclude that a more complicated galactic 
evolution model including radial flows, and possibly mass loss due to 
supernova explosions and winds, is necessary to explain the apparent 
chemical uniformity of the disk of NGC 1313. 
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1. INTRODUCTION 

NGC 1313 is often described as a galaxy in 
transition between the magellanic type galax- 
ies and normal disk galaxies such as the Sc 
galaxies M 33 (NGC 598) and NGC 300. 
Low- mass galaxies have maximum rotation 
velocity between 40 and 70 km/s, (Carig- 
nan & Freeman 1985; Carignan & Beaulieu 
1989; Cot, Freeman, & Carignan 1997), while 
normal galaxies of latest morphological types 
have values ranging from 90 to 110 km/s 
(Puche, Carignan, & Bosma 1990; Carignan 
& Puche 1990; Zaritsky et al. 1994). Chem- 
ically, Sc and Sd galaxies show very steep 
radial distributions of abundances (see e.g. 
Zaritsky et al. 1994), while low-mass irreg- 
ulars have very homogeneous abundance dis- 
tributions (Pagel et al. 1978; Kobulnicky 
k Skillman 1996; Roy et al. 1996; Devost, 
Roy, & Drissen 1997; Kobulnicky 1998). The 
galaxy NGC 1313 has a maximum rotation 
velocity very similar to that of NGC 300 and 
NGC 598. However, contrary to them, its 
radial distribution of oxygen abundance ap- 
pears very uniform, with a mean value around 
12 + log (O/H) = 8.4, a level similar to what 
is observed in the Large Magellanic Cloud 
(Pagel et al. 1980; Walsh & Roy 1997). Thus 
NGC 1313 is the most massive disk galaxy 
known showing very little or no radial gradi- 
ent in its radial abundance distribution. 

Low-mass galaxies are usually irregular in 
their optical appearance due to the smaller 
mass, and consequently, the smaller poten- 
tial prevents the formation of spiral structure. 
These galaxies have large relative quantities 
of atomic gas; because molecular gas is dif- 
ficult to detect in these systems, it is not 
yet understood if this is truely due to low 
content. The atomic gas distribution shows 
an exponential decrease with galactocentric 



radius from a central maximum. Moreover, 
many galaxies are undergoing a burst of star 
formation in the center (Greggio et al. 1992; 
Gallart et al. 1996; Skillman, Dohm-Palmer, 
& Kobulnicky 1998); they have low metallici- 
ties (Garnett et al. 1995). These facts and the 
intrinsic blue colors (Skillman 1996; Pildis, 
Schombert, & Eder 1997) indicate that low- 
mass galaxies are relatively unevolved. 

The rotation curve can generally be mea- 
sured far from the galactic center, and the ve- 
locity behavior suggests a large proportion of 
dark matter in these galaxies. However an al- 
ternate explanation within the framework of 
the MOND theory of gravitation, which gives 
flat rotation curves for the observable matter 
in galaxies, is possible (McGaugh & de Block 
1998). 

The systematic variation of abundances 
with galactocentric radius in normal spiral 
galaxies has been known for a long time (see 
references in Diaz 1989; Skillman 1997; Gar- 
nett et al. 1997; Henry 1998). Radial abun- 
dance gradients vary in amplitude. Corre- 
lations between their amplitudes and other 
galaxy characteristics have been clearly estab- 
lished (Vila-Costas & Edmunds 1992; Zarit- 
sky et al. 1994; Martin & Roy 1994). The 
radial gradient of abundances is steeper in 
later types of galaxies and the mean abun- 
dance level is moderate. For early type galax- 
ies, gradients appear flatter and mean abun- 
dances seem higher (Oey & Kennicutt 1993), 
but this is based on small samples due to the 
difficulty of measuring the O/H abundances 
in these systems. Gas distributions in normal 
galaxies also have some trends: the atomic 
gas shows an apparent hole in the central re- 
gion, while the maximum column density is 
reached further away into the disk. The po- 
sition of this maximum is shifted towards the 
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inner disk for late type galaxies, while for the 
early types, the maximum is further out in the 
disk (Molla, Ferrini, & Diaz 1997; Broeils & 
Rhee 1997). The proportion of gas in molecu- 
lar or atomic phases also depends on the Hub- 
ble type, with the larger ratios of molecular to 
atomic being the attribute of early type galax- 
ies (Young & Knezek 1989). We should point 
out that this relation has been obtained with- 
out taking into account the effects of chemical 
abundance and UV radiation field in the CO 
to H2 conversion factor. These effects have 
been demonstrated to be important in calcu- 
lating the molecular gas masses (see Verter 
& Hodge 1995; Wilson 1995). Molecular gas 
radial distributions have exponential shapes 
function of radius, mimicking the star 
formation rate distributions. Irregular galax- 
ies, with lower mass, luminosity, and absolute 
abundances, and exponential radial distribu- 
tions for gas and star formation rate follow 
the same behavior as above. Their only in- 
consistent feature with respect to those trends 
is their flat radial abundance distributions. 

2. CHEMICAL EVOLUTION MOD- 
ELS OF DISK GALAXIES 

The multiphase chemical evolution model, 
developed by Ferrini et al. (1992; 1994), has 
been applied to several spiral galaxies (Molla, 
Ferrini & Diaz, 1996: hereafter MFD96) with 
different morphological types. The model re- 
produces the observed correlations by varying 
the characteristics infall rate and the cloud 
and star formation efficiencies with the galaxy 
morphological type and/or the arm class as 
defined by Elmegreen & Elmegreen (1987). 
The observed radial distributions of atomic 
and molecular gas, of oxygen abundances and 
of the star formation rate are used as ob- 
servational constraints. We can infer from 



the previous simulations that a galaxy less 
massive than NGC 300 and NGC 598, which 
have been modeled by MFD96, would have 
evolve more slowly, because of lower efficien- 
cies in forming molecular clouds and stars; 
this would produce an exponential radial dis- 
tribution of atomic gas and a steeper abun- 
dance radial gradient than those obtained in 
MFD96 for NGC 300 and NGC 598. The 
strong gradient would be due to the high star 
formation rate in the central regions. Mas- 
sive stars from this episode would eject a 
elements, increasing the oxygen abundance 
mainly in this zone, while the rest of the disk 
would maintain a low value. The difference 
between the central level of abundance and 
the value in the outer region would give the 
appearance of a steep radial gradient. 

An important parameter influencing the 
shape of the radial distribution of abundances 
is the initial radial mass distribution, or in- 
directly the total mass of the galaxy. The 
radial mass profile, M(R), may be indeed dif- 
ferent for galaxies with the same Hubble type, 
allowing different abundance distributions or 
gradients. Our aim, in this work, was to ver- 
ify whether the multiphase chemical evolu- 
tion model would actually predict a steep ra- 
dial gradient for the oxygen abundance in the 
galaxy NGC 1313 or could produce the ob- 
served flat gradient using a range of different 
assumptions, while reproducing the observed 
radial distributions the gas and star forma- 
tion. 

Not yet included at present in our model is 
the presence of a bar. NGC 1313 has a short, 
but rather strong bar (Martin 1995). Barred 
galaxies have a flatter radial gradient (Ed- 
munds & Roy 1993; Martin 1992; Zaritsky, 
Kennicutt, & Huchra 1994; Martin & Roy 
1994, 1995; Ryder 1995; Roy 1996). It has 
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been shown that bars induce large-scale mix- 
ing due to radial flows (Friedli, Benz, & Ken- 
nicutt 1994; Friedli & Benz 1995). The radial 
flows produce a net accumulation of mass in 
the inner regions of the disk, and may enhance 
a central star formation. This results first in a 
flatter abundance gradient for the outer disk 
and in a steepening of the gradient for the in- 
ner disk; this behavior is a local and transient 
phenomenon which lasts a few hundred mil- 
lions years. On the longer timescale, radial 
mixing dominates and the global gradient is 
made flatter than before the bar formation. 
The dynamical and morphological evolution 
of low mass galaxies has not been studied yet 
by theoreticians because of the difficulties in 
building realistic codes of disk galaxies where 
the gas component dominates the mass (D. 
Friedli, private communication). 

The multiphase model uses as input the 
total mass radial distribution, obtained from 
the rotation curve; the possibility of the rota- 
tion curve changing with time or the presence 
of dynamical processes, such as radial flows of 
gas, are not taken into account. The model 
has been already applied to three massive 
barred galaxies (Molla, Hardy & Beauchamp 
1998), and we know that it is not possible to 
fit the abundance data in the bar region with 
the multiphase approach. To describe these 
zones, a chemical evolution model including 
radial flows of gas must be developed. The 
galaxy NGC 1313 has a bar, but it will be 
assumed to be small enough not to have an 
important effect. This paper assesses the va- 
lidity of this simple assumption. 

3. PROPERTIES OF NGC 1313 FROM 
THE MULTIPHASE MODEL 

We have applied the multiphase model to 
the disk galaxy NGC 1313 which is an iso- 



lated very late-type barred galaxy. Its mor- 
phological type lies between irregular magel- 
lanic spiral (Sm) and the latest type of normal 
spirals (Sd). We modeled the galaxy chemi- 
cal evolution as in MFD96. We assumed that 
the protogalaxy is a spheroid of primordial 
gas. The total mass included in the spheroid 
and its radial distribution are inferred from 
the rotation curve. The rotation curve is de- 
rived from the H I 21 cm line data of Ryder 
et al. (1995); the maximum rotation velocity 
given by these authors is V max ~ 110 km/s. 
The rotation curve shows an apparent hump 
at 7 kpc, but the authors warn that it may be 
an artifact due to the high inclination of the 
galaxy. We use their fit to the data - with a 
smoother behavior - to derive the total mass 
of the corresponding spherical protogalaxy. 

The galaxy is divided into 10 concentric 
cylindrical regions 1 kpc wide; a .halo zone 
and a disk zone correspond to each galac- 
tocentric step. As the gas begins to form 
stars in the halo, it collapses onto the equa- 
torial plane, forming the disk. The collapse 
proceeds at a rate that is inversely propor- 
tional to the collapse time scale. Due to 
its dependence on the total local mass sur- 
face density, this timescale, r(R), changes as 
r oc ct(R)- 1 / 2 (Jones & Wyse, 1983). We as- 
sume this radial dependence to have an expo- 
nential form: t(R) = roe -11 / 1 , where / is the 
disk scale length taken as 2.5 kpc (Ryder et al. 
1995). The characteristic collapse timescale, 
T , depends on the total mass of the galaxy 
through the relation t oc M 9 T (Gallagher, 
Hunter & Tutukov, 1984), where M 9 is the to- 
tal mass of the galaxy in 10 9 M , and T is 
its age, assumed to be 13 Gyr. We calculate 
the collapse timescale r with the expression 
t = T (M 9i g a i/M 9i MWG)~ 1/2 , where the sub- 
script MWG refers to the Milky Way, and r 
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Table 1: 


NGC 1313 Model Input Parameters 




Name 


Meaning 


Value 


To 


Collapse Time Scale 


12 Gyr 


ex 


Efficiency for the star formation in the halo 


0.009 




Efficiency for the cloud formation in the disk 


0.05 




Efficiency for the spontaneous star formation in the disk 


0.003 




Efficiency for the stimulated star formation in the disk 


2.5xl0~ 8 



taken as 4 Gyr, is the collapse timescale used 
for the solar neighborhood region (Ferrini et 
al. 1994). Total masses in units of M 9igal 
are calculated as M gal = 2.32 x 10 5 R ga iV^ ax 
(Lequeux, 1983), with R ga i = 15 kpc; V max 
has the value 110 km/s already quoted. The 
collapse time finally selected is r co u = 12 Gyr; 
it is normalized to a region equivalent to the 
solar region, and is located at a galactocen- 
tric distance Ro which is at 2 kpc based on 
a value of 1.73 kpc for the effective radius of 
NGC 1313 (Ryder et al. 1995). 

The initial mass function (IMF) is that of 
Ferrini, Palla, & Penco (1990), and the nu- 
cleosynthesis prescriptions are the same as 
used in other multiphase models: Woosley 
& Weaver (1986) for massive stars; Renzini 
& Voli (1981) for low and intermediate mass 
stars and Nomoto, Thielemann & Yokoi (1984) 
and Branch & Nomoto (1986) for type I su- 
pernova explosion yields. 14 N production is 
calculated with Maeder (1983) prescriptions. 
The star formation in the halo is assumed to 
be a Schmidt law proportional to the gas den- 
sity (smooth distribution of atomic gas) with 
n = 1.5 as an exponent. The proportional- 
ity coefficient has a non linear dependence on 
the radius through the halo volume, for ev- 
ery zone or region, and on the efficiency of 
the process, ex, which is assumed constant 
for every region and galaxy. 

Star formation in the disk is considered as 



a two step process: first, molecular clouds 
form from the uniform atomic gas; then stars 
are created by cloud-cloud collisions (sponta- 
neous star formation). The induced star for- 
mation from the interaction of massive stars 
(through their winds and/or supernova shells) 
with molecular clouds is also included (stimu- 
lated star formation). The observed star for- 
mation in molecular clouds provides the em- 
pirical basis for our star formation law, which 
may be closer to the true physical mechanism 
than other assumed laws. The star formation 
rate depends on the masses of the different 
phases of the matter involved with propor- 
tionality coefficients that depend on some ef- 
ficiency or probability factors (Table 1). 

The induced star formation is proportional 
to the efficiency, e a , which is usually assumed 
constant for all regions and galaxies. The 
two other processes of formation have 
dial dependence due to the volume effect, and 
a dependence on the efficiencies e M for the 
cloud formation and e# for the stimulated 
star formation. The range for their charac- 
teristic values is selected by taking into ac- 
count the Hubble type, and/or the arm class 
for the galaxy and the previous values for all 
the Hubble types (see Table 4 of MFD96). 
For NGC 1313, these values are e M = 0.05 and 
€h = 0.003 and are applied to the equivalent 
solar region with radius Ro- 

Final values for the efficiencies and collapse 



5 



0.5 



X 

CD 

ao 
o 

+ 



7.5 



1 i 

Model 1 
Model 2 
Model 3 
Model 4 




4 6 
Radius (kpc) 



1 



Fig. 1. — Predicted radial distribution of the oxygen abundance 12 + log (O/H) at the present 
time for NGC 1313: Model (1) is represented by the short dashed line, model (2) by the solid line, 
model (3) by the long dashed line, and model (4) by the dotted line. Observational data points 
are from Hn region O/H abundance measurements by Walsh & Roy (1997). 



times have been selected to construct a series 
of models fitting the observed radial distribu- 
tions of atomic gas (Hi), star formation rate 
surface densities and oxygen abundance of the 
gas phase. Observational data for H I surface 
densities are from Peters et al. (1994) and 
Ryder et al. (1995). The star formation ra- 
dial distribution is derived from the Ha flux 
radial profile of NGC 1313 by Ryder (1993); 
we normalized it to the equivalent to solar re- 
gion value ^/ . Oxygen abundances are those 
of Walsh & Roy (1997), who derived O/H 
abundances using the calibration of the neb- 
ular diagnostic line ratio ([On] + [Olll]/H/3 
by Edmunds & Pagel (1984). Unfortunately, 
there are no data on the molecular gas. 

We ran several models trying to reproduce 
the observed radial distributions in a self- 
consistent manner. We show the results for 
four different models in figures 1 to 4. These 



different models have almost the same effi- 
ciencies e M and e# at the distance equiva- 
lent to solar radius (2 kpc), but use differ- 
ent assumptions for radial dependences, as we 
will explain in the next paragraph. Figure 1 
show the predicted radial abundance distri- 
bution compared with the corresponding ob- 
servational data; the black dots are the Hn 
region O/H measurements of Walsh & Roy 
(1997). Figure 2 displays the theoretical ra- 
dial distributions of the star formation rate 
surface density, normalized to the value ^q; 
the actual rates, from the azimuthally aver- 
aged Ha flux, are shown by the black dots. In 
Figures 3 and 4, the inferred radial distribu- 
tions of atomic H I gas and molecular H 2 sur- 
face densities are shown along the true data as 
squares or black dots for H I from Peters et al. 
(1994) and Ryder et al. (1995), respectively. 
Model 1 is the standard multiphase model. 
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Fig. 2. — Predicted radial distribution of the star formation rate surface density, normalized to 
the value \l/o- Observational data are from Ha flux measurements by Ryder (1993). The lines 
referring to different models have the same meaning as in Figure 1. 



The efficiency to form molecular clouds e M 
varies with galactocentric radius, being high- 
est for the inner regions of the disk. The effi- 
ciency 6h to form stars from cloud-cloud col- 
lisions process is kept constant with the ra- 
dius. This model (1) is shown in Figures 1 to 
4 with the dashed line. The atomic gas distri- 
bution reaches a maximum at 4 kpc (Figure 
3); the O/H gradient is very steep (~ —0.128 
dex/kpc) (Figure 1). This is an expected re- 
sult in light of the sequence of models pre- 
viously run for several other galaxies. Actu- 
ally, this model is very similar to the models 
of MFD96 for NGC 300 and NGC 598 using 
the same multiphase approach; NGC 300 and 
NGC 598 are observed to have maxima of H I 
in the inner disk (at ~ 3 kpc ) and steep radial 
O/H gradients. However, NGC 1313 is not 
observed to have these features. We need a 
model giving an exponential atomic gas dis- 
tribution, and at the same time, a flat ra- 



dial abundance gradient. Though successful 
for two galaxies of similar mass and luminos- 
ity, our calculations fail to adjust the relevant 
physical distributions in NGC 1313. 

Our assumptions for Model 1 do not ap- 
pear valid. We have supposed that the effi- 
ciency to form molecular cloud, e M , is larger 
in the inner disk due to the effect of the spi- 
ral density wave. Taking into account that 
NGC 1313 is more irregular in its optical ap- 
pearance than NGC 300 and NGC 598, the 
two galaxies of comparison, we may relax this 
assumption and check whether a model with 
an efficiency constant along the radius gives 
a better result. 

We also assumed that once clouds are formed, 
the efficiency to form stars from molecular is 
constant. But, with this assumption, the ob- 
served exponential function of the star for- 
mation rate was not fitted for NGC 300 nor 
NGC 598, because the distribution had the 
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tendency to flatten in the inner disk. To al- 
leviate this problem, recent models (Molla, 
Hardy & Beauchamp 1998) have been com- 
puted with smaller efficiencies at shorter radii. 
The reason for this variation is that gas den- 
sities need to be above a threshold density for 
star formation to proceed (Kennicutt 1989). 
This critical density depends on the disper- 
sion velocity, usually constant along the disk 
in spiral galaxies, and on the epicyclic fre- 
quency. Although the velocity dispersion usu- 
ally does not vary much across the disk, it 
is possible to find disk regions where it is 
larger. In normal galaxies, these zones may be 
those where the bar, or the bulge, introduces a 
strong perturbation in the velocity field. Ry- 
der et al. (1995) have measured dispersion ve- 
locities in NGC 1313, where they average 14 
km/s; they are larger than in normal galaxies. 
They also found that the velocity dispersion 
reaches a value of 25-30 km/s near the center 
of NGC 1313. If star formation efficiencies is 
higher where the dispersion velocity is lower, 
stars should form relatively more easily in the 
outer disk. This would be consistent with our 
hypothesis that it is more difficult to form 
stars in the center of NGC 1313 than in the 
comparison galaxies. 

Both these assumptions has been used to 
construct Model 2, shown by solid lines in 
Figures 1 to 4. Due to the larger efficien- 
cies to form stars further away in the disk, 
we are successful at getting a flat O/H radial 
distribution. The radial distribution of the 
atomic gas is also satisfactorily reproduced; 
the molecular gas forms in the same propor- 
tion at all radii, by exhausting the atomic gas 
in the outer regions of the disk, where the to- 
tal disk mass is already small. Where Model 
2 fails is in the large star formation rate pre- 
dicted in the outer radii in comparison with 



the Model 1. A flat radial gradient of abun- 
dances is obtained, but the model does not 
reproduce the rate inferred from the observed 
Ha profile (Figure 2): the predicted star for- 
mation rate is flatter than observed, with un- 
realistically high values at larger radii. 

Both these assumptions has been used to 
construct Model 2, shown by solid lines in 
Figures 1 to 4. Due to the larger efficien- 
cies to form stars further away in the disk, 
we are successful at getting a flat O/H radial 
distribution. The radial distribution of the 
atomic gas is also satisfactorily reproduced; 
the molecular gas forms in the same propor- 
tion at all radii, by exhausting the atomic gas 
in the outer regions of the disk, where the to- 
tal disk mass is already small. Where Model 
2 fails is in the large star formation rate pre- 
dicted in the outer radii in comparison with 
the Model 1. A flat radial gradient of abun- 
dances is obtained, but the model does not 
reproduce the rate inferred from the observed 
Ha profile (Figure 2): the predicted star for- 
mation rate is flatter than observed, with un- 
realistically high values at larger radii. 

In Models 3 and 4, variation in the rela- 
tive importance of stimulated star formation 
was tried. This is done by changing the pa- 
rameter e a which measures the rate of interac- 
tion between massive stars and the surround- 
ing molecular clouds. Normally, this param- 
eter is assumed constant, because this it is 
presumed to be a local process (the star for- 
mation rate depends on the local gas surface 
density) restricted to a relatively small re- 
gion. However, this process probably depends 
on the ultraviolet luminosity and on the den- 
sity of the gas and dust that may shield the 
region from the radiation field. To explore 
the effect of enhanced stimulated star forma- 
tion, we decreased the efficiency of sponta- 
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Fig. 3. — Predicted radial distribution of the atomic gas Hi surface density. The observed data 
are taken from Ryder et al. (1995) - circles - and Peters et al. (1994) - squares. The lines 
referring to different models have the same meaning as in Figure 1. 



neous star formation to let stimulated star 
formation dominate. If this quantity has a 
radial distribution with an exponential func- 
tion, the star formation rate can be adjusted 
with this kind of model. Model 3, having a 
low efficiency e a , leads to an exponential be- 
havior for the star formation rate and a flat 
abundance gradient, but the level of abun- 
dance, 12 + log(0/H) ~ 7.7, is far below the 
observed level. If we increase the star for- 
mation rate to raise the abundance level by 
increasing the efficiency e a or the efficiency 
- this is Model 4 -, the star formation rate be- 
comes again unrealistically large in the outer 
parts. 

4. DISCUSSION AND CONCLUSIONS 

We tried to reproduce the radial abun- 
dance distribution observed in the ionized in- 
terstellar gas of the galaxy NGC 1313, which 



has a flat abundance gradient; we compared 
with models for two similar disk galaxies, 
NGC 300 and NGC 598, showing a very steep 
radial distribution of atomic gas. We failed 
for NGC 1313, because we were unable to fit 
the present observed star formation rate ra- 
dial distribution in a self-consistent way with 
other related distributions, those of gas and of 
O/H. Taking into account the observed global 
radial abundance gradients, the radial distri- 
butions of star formation rates and of gas sur- 
face densities, we predicted very steep radial 
O/H distributions for all three galaxies; the 
observed distributions for atomic gas and star 
formation rate follow our expectations, but 
that of O/H does not in NGC 1313. 

Edmunds & Roy (1993) remarked that global 
O /H radial gradients seem to disappear when 
the spiral structure no longer exists. Explor- 
ing further this idea, we ran a model with a 
constant and low value for the efficiency of 
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Fig. 4. — Predicted radial distribution of the molecular gas H 2 surface density. The model lines 
have the same meaning as in Figure 1. 



molecular cloud formation and a small value 
(decreasing in the inner regions) for the effi- 
ciency of spontaneous star formation, to try 
to get a flat gradient in the radial abundance 
distribution. It did not work: a flat gradient 
was obtained, but the predicted star forma- 
tion radial profile was much flatter than ob- 
served. 

On the other hand, the star formation rate 
being highest at the center, we expect the 
oxygen abundance to be maximum there, as 
a consequence of enhanced production of ele- 
ments. Could then the effective enrichment 
be less efficient because of loss of metals, 
thus resulting in a flatter gradient? Various 
possibilities for element loss exist. Tenorio- 
Tagle (1996) has already suggested that re- 
cent heavy element enrichment is unobserv- 
able because the fresh ejecta are hiding in 
a hot phase. A second possibility involves 
mass loss due to collective winds produced 
by supernova explosions (Charlton & Salpeter 



1989). The mass ejected expel a-elements 
produced in the first stellar generations out 
of the disk. This hot gas condenses later and 
falls back over a long time-scale, resulting in 
an homogeneous oxygen abundance distribu- 
tion in the interstellar gas. Actually, both 
possibilities may combine: supernova explo- 
sions produce a galactic wind which ejects a 
proportion of elements already created. The 
hot gas diffuses away from its injection loca- 
tion and eventually falls on any region of the 
disk. It eventually mixes with the cold phase, 
and the elemental abundances seen in H II re- 
gions show a uniform radial distribution. This 
possibility is consistent with the calculations 
of Hunter, Elmegreen & Baker (1998) in their 
study of star formation mechanisms in irreg- 
ular galaxies. The authors suggest, as a cri- 
terion for the onset of star formation, a criti- 
cal density allowing the existence of cool gas; 
this critical density H t h depends on the to- 
tal pressure and on the average stellar radi- 
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at ion field. They found that most galaxies 
have E g /Et/j < 1, which suggests that the gas 
is in the warm phase in irregular galaxies due 
to the lower pressure. It means that it is hard 
to sustain a cool phase in the Hi phase, and 
that most of the gas is warm. 

A big caveat is in order: why would the loss 
mechanisms work differently in NGC 1313, 
as compared to NGC 598 or NGC 300 which 
have about the same total mass? Could this 
be due to a different merger history? The 
supershell found by Ryder et al. (1995) in 
NGC 1313 might be explained as the effect 
of a merger. However Ryder et al. demon- 
strate that the shell must originate from an 
internal event. Furthermore there is no indi- 
cation from the H I velocity field nor from the 
photometry of the disk that a merger event 
took place in the recent history of NGC 1313. 
Thus the possibility of a merger as the cause 
of the flat gradient is negligible. 

Furthermore, flat gradients or uniform abun- 
dance distributions are a normal feature of 
low-mass irregular galaxies, not an excep- 
tional characteristics. We must find a gen- 
eral mechanism producing a flat oxygen abun- 
dance distribution, and giving at the same 
time consistent corresponding radial distribu- 
tions for star formation and atomic gas sur- 
face densities. This mechanism should also 
work for other low mass irregular galaxies 
such as those in Table 3 of Walsh & Roy 
(1997). From a statistical point of view, 
the action of radial flows induced by bars 
as the cause of flat gradients is more prob- 
able than rare merger events. Radial flows 
of gas appear as a consequence of the non- 
axisymmetric gravitational potential; these 
radial flows are efficient at mixing the inter- 
stellar medium and can produce a flat gra- 
dient over a timescale of ~1 Gyr (Friedli, 



Benz, & Kennicutt 1994; Edmunds & Green- 
how 1995). 

In the case of the galaxy NGC 1313, sev- 
eral authors (Blackman, 1981; Marcelin & 
Athanasoula 1982; Ryder et al. 1995) came 
to the conclusion that the velocity field pat- 
tern could only be explained by introducing 
a radial velocity component. This is typi- 
cally the case for barred galaxies; the radial 
component may be as small as 20 km/s or 
as large as 80 km/s. For galaxies with big- 
ger bars than NGC 1313, radial flow veloc- 
ities are ~ 40 km/s (Blackman 1981; Ball 
1992). For the Large Magellanic Cloud which 
is smaller than NGC 1313, the recent work of 
Kim et al. (1998) reveals streaming motion 
in the vicinity of the stellar bar; from these 
data, we infer the amplitude of the pertur- 
bation in the velocity field to be ~ 10 — 15 
km/s. Exact calculation of the velocity of ra- 
dial flows from the observed H I velocity field 
requires a complex numerical simulation code; 
however by scaling for galaxy mass and bar 
strength, we estimate that the radial compo- 
nent due to a bar may be 15 — 20 km/s for 
NGC 1313. Assuming 20 km/s, chemical mix- 
ing of the gas may be achieved in a timescale 
between 0.5 and 2 Gyr following the analysis 
of Martin & Roy (1995) and Roy & Kunth 
(1995). The shortest timescale is obtained by 
simply dividing the radius of the disk by the 
mean radial flow velocity; however this is too 
simplistic. The longest timescale is obtained 
by calculating the pseudo-diffusion timescale 
t = R 2 / (v I), where R is the radius of the disk, 
v the mean radial flow, and I is the width of 
the annuluses over which the net radial flow 
is either inward or outward (see Fig. 2 in 
Friedli, Benz, & Kennicutt 1994). This up- 
per value is probably the most plausible, and 
is much shorter than the Hubble time. Thus 
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radial flow should influence the radial abun- 
dance distribution if the galaxy is several Gyr 
old. Also with radial flows, a large quantity 
of gas would accumulate in the center of the 
galaxy, and could generate a starburst. This 
is consistent with soft X- ray emission from 
the center of NGC 1313, which could be pro- 
duced by a low-luminosity starburst or AGN 
(Colbert et al. 1995). 

Because of mass conservation, radial flows 
must be compensated by the formation of 
stars in the center which acts sink. If 
not, the disk will deviate from a 1/r profile 
and will not be stable (Struck-Marcell 1991). 
Massive star formation will insure that the 
excess of mass is expelled out of the region 
through some sort of wind phenomenon sim- 
ilar to that produced by supernovae explo- 
sions. Maintaining a 1/r profile imposes a 
limit on the radial flow velocity; this limits 
depends on the total local surfacic mass den- 
sity and on the star formation rate in the cen- 
tral region (Struck-Marcell 1991). We have 
calculated this limit to be v r = 0.2 km/s for 
NGC 1313. If the bar induces a radial flow 
with a velocity larger than ~1 km/s, there 
should be some mass loss. Therefore, both 
radial flows and mass loss by winds must be 
incorporated into chemical evolution models. 
We also need a multi-zone model to calcu- 
late radial dependences and exchanges be- 
tween radial zones. We hope to show in a near 
future results from models with these charac- 
teristics. 
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